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In a recent article- Yuen-Zhuo et al. (YZ) have 
proposed an extension of our theorem of Stochastic 
Time-Dependent Current Density Functional Theory 
(STDCDFT)^ The main claim of this generalization 
is a mapping between open and closed quantum systems, 
namely that the dynamics of the current and particle 
densities of an open quantum many-body system can be 
uniquely determined by the time evolution of the same 
quantities of an effective non-interacting closed system. 
If this were true, there would be no need to consider the 
evolution of a density matrix via a master equation, but 
rather the dynamics of the current and particle densities 
could be obtained via the time evolution of the standard 
Schrodinger equation. From a computational point of 
view this would be dramatically cheaper than solving for 
a master equation^ and would also leave open the possi- 
bility of analytical solutions, asymptotic expansions, etc. 
In this comment we point out some serious deficiencies 
in the proof of YZ's theorem that invalidate their main 
results. 

Let us begin with a critical examination of Theorem 1 
in Ref. To do this, let us recall that in the theory of 
open quantum systems, the equation of motion for the 
ensemble-averaged particle density, n(r,t), is given by- 



d t n(r,t) = -V -j(r,t)+T B (r,t) 



(1) 



where Tb ( r , t) describes the density modulation induced 
by the presence of the bath and j(r,t) is the ensemble- 
averaged current density^ Eqn 1 can be obtained from 
the master equation for the density matrix, p 



dt P {t) 



H(t), P (t) 



dt'K(t,t')p(t') 



(2) 



Here, K(t,t') is a memory kernel that describes the ac- 
tion of the bath on the system, and H(t) is the Hamil- 
tonian of the system that evolves under the action of 
an external vector potential possibly time dependent (we 
work, as in YZ's work in a gauge in which the scalar 
potential has been set to zero). 

Theorem 1 in Ref. [l| then states that there is a one- 
to-one correspondence between a vector potential and 
the pair of functions given by the current and particle 



densities. However, this statement cannot be correct 
as a simple mathematical characterization of the func- 
tional spaces connected by this mapping shows. If, as 
YZ seem to argue, the particle and current densities are 
independent functions, then they define a 4-dimcnsional 
functional space. This amounts to saying that Tb{t, t) 
is not completely determined by the sole knowledge of 
n(r,t) and j(r,t). Since a general vector potential is 
a vector of three independent functions, Theorem 1 in 
YZ's work would imply that a 4-dimensional functional 
space formed by the density plus the three components 
of the current density is locally homeomorphic to a 3- 
dimensional functional space spanned by the vector po- 
tential. This is obviously incorrect. A similar problem 
is present in the connection between the scalar potential 
and the current density^ The solution to this inconsis- 
tency is that the continuity eqn (1) is indeed a non- linear 
equation in the particle and current densities from which 
the particle density can be obtained starting from the 
current density, once the bath operator and initial con- 
ditions have been fixed^^ This amounts to saying that 
3~B{ r , t) is a functional of n(r, t) and j(r,t), or better 
of j(r,t) alone, and that equation (1) admits a unique 
physical solution. 

It is also worth pointing out that our proof of the 
theorem^ as well as that of the theorem of standard TD- 
CDFT^ does not require that the particle density in the 
auxiliary system is equal to the particle density in the 
original one£ In those proofs one only needs to determine 
the n-th time derivative of all quantities, and only the 
(n+l)-th time derivative of the vector potential. The n- 
th time derivative of the particle density is obtained from 
eqn (1) (or, in a closed system from the same equation 
with Tb = 0). That equation, however, does not contain 
any (n+l)-th time derivative. Therefore, the equation 
for the vector potential in the auxiliary system we use 
in our proo£ is still a recursive relation, with a unique 
solution provided the initial conditions. The theorem 
then guarantees the one-to-one correspondence between 
external vector potential and ensemble-averaged current 
density, leaving quite open the (possibly difficult) task 
of obtaining the ensemble-averaged density from the cur- 
rent density, when the continuity equation is given by 
eqn (1). 
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Apparently, the motivation of YZ's work to extend 
our STDCDFT theorem appears connected to this point, 
namely that eqn (1) may not always be satisfied, i.e., that 
the contribution Tb (?", t) is not equal in the real and KS 
systems, if one keeps the bath operator fixed [see discus- 
sion after eqn (18) in their manuscript]. Besides the fact 
that this statement is not mathematically proven - and 
following the discussion of the previous paragraph, not 
strictly necessary - we could argue that, since in calculat- 
ing the particle and current densities one fixes only a few 
degrees of freedom of the density matrix, two density ma- 
trices can indeed produce the same average particle and 
current densities (in two different systems) and satisfy 
eqn (1) with the same bath operator. In addition, un- 
like what YZ claim, we have never argued that the term 
Tb (r, t) needs to be smalts or to vanish identically. In- 
deed, when discussing the continuity equation (1) in our 
formalism (sec. II-D of Ref. 0) we have requested that it 
uniquely fixes the ensemble-average particle density once 
the average current density is given, with or without the 
contribution Tb { r i t) ■ This is definitely true for the initial 
conditions, when the system and bath are (by hypothe- 
sis) uncorrelated and thus Tb(t, to) = 0. In our original 
paper we have then separated the following discussion in 
two cases. On the one hand we have considered that this 
extra term identically vanishes at any instant of time. 
This happens, for example, when we consider local bath 
operators. In this case we do not need to go any fur- 
ther because the density can be obtained uniquely from 
the current density in the standard way. On the other 
hand, if the extra term Tb (r, i) does not vanish identi- 
cally, the only physical option left is to assume that eqn 
(1) uniquely determines n(r, t) given j(r, t), once the ini- 
tial conditions and bath operator are fixed. Indeed, in 
our proof we construct the average density and current 
density in time, once the initial conditions are assigned, 
and show that the term Tb (r, t) is some functional of the 
current and particle densities. At this point, the (now 
possibly nonlinear) equation of motion (1) for the av- 
erage particle density admits one and only one physical 
solution given the average current density, once the initial 
state has been assigned^ 

We now want to critically examine Theorem 3 in Ref. [l|, 
where the main result of that paper is presented: YZ 
claim to show that it is possible to construct a unique 
closed non-interacting quantum system that mimics the 
dynamics of n(r,t) and j(r,t) of the real open interact- 
ing system. We want to point out that, following YZ's 
reasoning, one can in fact find many (possibly infinite) 
closed non-interacting quantum systems that reproduce 
the dynamics of the exact current and particle densities. 
However, this is in contrast with Theorem 1, that claims 
that given the densities only one such system should ex- 
ist. In their proof, at step 1, YZ state that the vector 
C(r,t) can be uniquely determined by their eqn (23) 



once a boundary condition in space that fixes an arbitrary 
function of time is assigned. However, the above equation 
is not the unique solution to eqn (22) of Ref. 0, which 
reads 

dn(r,t) ~ -.. , ( eC(r,t) . A 

— i^i = _V-j(r,i)-V- I — i-Li n (r,t)J . (4) 

For instance, easily fulfilling the assigned boundary con- 
dition, one can add the curl of an arbitrary vector, g(r, t), 
to the function C(r,t) and still satisfy their eqn (22). 
Indeed it is easily proven that C'(r,t) = C(r,t) + (V x 
g(r, t))/n(r, t) still satisfies eqn (22) if C(r, t) satisfies the 
same equation. If we now continue with the proof and use 
C'(r, t) we arrive at a new vector potential A' KS (r, t) that 
gives the same n{r,t) and j(r,t) as the couple Aks(t, t) 
and C(r, t). This, however, is in contradiction with Theo- 
rem 1 that claims a one-to-one mapping between n(r,t), 
j{r,t) and Aks{t, t), when all the other operators and 
the boundary condition, are kept fixed." This ambiguity 
reflects the fact that YZ are trying to mimic the effect of 
a scalar function - the term TB(r,t) in eqn (1) - with a 
vector function, C(r,t), without imposing strict bound- 
ary conditions (BCs) on C(r,t). 

A simple solution to the aforementioned problem may 
appear by setting VxC(r,f) = 0. However, it is im- 
portant to realize that the imposition of certain bound- 
ary conditions on the dynamics of these quantities has a 
direct impact on the uniqueness of the results. For ex- 
ample, assuming that C(r,t) reaches a certain uniform 
limit when \r\ — > oo, might be inconsistent with the bath 
operator acting on the true many-body system. Indeed, 
certain bath operators can be strongly non-local in space, 
effectively transferring charge from one region of space 
to another, with the two arbitrarily far from each other. 
Therefore, the BCs on V x C(r, t) have not a clear phys- 
ical origin or relation to any physical observable. Fix- 
ing their value to obtain one solution appears utterly 
arbitrary. Notice that a similar problem appears also in 
the case of the standard theorem of TDDFT (see for ex- 
ample Ref. 13j) where the additional boundary condition 
n(r)V AV(r,t) —> when \r\ — > oo is added to the proof. 
However, in this case, while in principle this condition 
is arbitrary and one may choose another condition, this 
choice is motivated by physical arguments that are valid 
for a wide range of systems. The same considerations 
instead do not apply to all the components of the vector 
C(r,t). Therefore, the proof of the theorem as it is for- 
mulated in Ref. [l| cannot hold for general bath operators. 

Finally, we want to comment on a fundamental but 
important issue. YZ's initial assumption for their theo- 
rem is a closed equation of motion for the density ma- 
trix. As we have discussed at length in our previous 
publications^^ this is not a solid starting point for a for- 
mulation of DFT for open quantum systems. This is due 
to both the possible loss of positivity of the density ma- 
trix if an equation of motion of such quantity is employed 
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with the Hamiltonian and/or bath operator (s) dependent 
on time^ and the fact that the KS Hamiltonian does de- 
pend on internal degrees of freedom. Starting from the 
master equation formulation of the same problem, one 
needs to exclude from the outset the possibility that the 
Hamiltonian of any auxiliary system with different in- 
teraction potential (and hence the KS Hamiltonian) de- 
pends on the internal degrees of freedom. Otherwise, for 
such a system no closed density-matrix equation can be 



obtained. In other words, one needs to start from an hy- 
pothesis that constitutes part of the final thesis. It is only 
when one starts from a stochastic Schrodinger equation 
for the state vector that one can prove that the exact KS 
Hamiltonian depends only on the average current density. 
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